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Abstract 



By examining both the divergence of the velocity vector in orthogonal Carte- 
sian coordinate space T of dimension W and the structure of the Hamil- 
tonian determining a system trajectory, it is shown that the standard Liou- 
-. I ville equation cannot describe anything more than linear motion for typical 

^^ ■ Hamiltonians with separated momentum and space variables and some sig- 

0^ , nificant consequences such as the Poincare recurrence theorem do not obtain 

for such Hamiltonians. A new stochastic equation which is everywhere in 
^3 ' principle discontinuous is developed for dynamical systems described by a 

general Hamiltonian which is a functional of the space and momentum vari- 
ables, and where the average trajectory of a system point is proven to be 
fi ■ orthogonal to any constant energy surface consonant with the system energy 

^ . at equilibrium. This equation does not assume the presence of binary colli- 

>v>( I sion only, as required in the standard first-order Boltzmann equation, and is 

^ ' therefore suitable to describe dense systems as well, and may be viewed as 

an alternative to the latter. Some new macroscopic variational principles for 
non-equilibrium thermodynamical systems are proposed, where one object for 
future work would be to relate the microscopic description given here with 
the macroscopic principles. 

1 Introduction 

The classical Liouville equation of statistical mechanics is described by trajec- 
tories in so-called T space. In this 7 (or F) space, N* systems are represented, 
each with coordinates (p,q), where the system has N particles, each of which 
has f degrees of freedom for the momentum and space variables. This leads to 
a 2fN dimensional T space where each system may be represented by a point 



traveling with a velocity v = (p,q). From the definition of the density-in- 
space p given by p = jyt , , the following consequences obtain: 



dp 
where 



^^+dlv{p^,)=Q (1) 



The fortuitous cancelation of terms ||l| in the 2nd R.H.S. term of Eq.(2) 
from the relation qj = ^—,Pj = ~3~ involving the Cartesian orthogonal 
coordinates lead to div v = or 

t = -E (^^-^ +^^-^) = -"•^^"'^^ (3) 

Eq. (3) is the Fundamental Liouville equation. Because of the fortuitous can- 
celation of terms involving the momentum and position coordinates, it is not 
immediately apparent that Vi = constant for each component Vi of v and 
in fact it is assumed in general as a fundamental fact that v is not a con- 
stant vector, and so this assumption must be tested by testing the divergence 
independently. Furthermore, the density p{p, q) is sometimes written inde- 
pendently of the positional coordinates, as in the case of a perfect fluid inside 
a region with a constant potential, such as in the Gibbs' Canonical ensemble 
written p{p) = Cexp —(3H{p, q) where it cannot be assumed that a fortuitous 
cancelation of terms occurs because of the incomplete variables specified, so 
another method must be utilized to check for the velocity values. Because the 
Hamiltonian is assumed not to be explicitly a function of the time t then the 
divergence would involve taking limits into regions oi p — q space which are 
not constant in energy, whereas a particular system trajectory would always 
conserve energy. Hence the macroscopic behavior of the system predicted 
from the density is dependent ultimately on the nature of the density flow in 
r space, and hence on the divergence of the vector field. 



2 Test of divergence term and its conse- 
quences 

The divergence is computed from the traditional definition 



div3 = lim / J.dS /5V (4) 




where 5V is the volume increment Jqr^fc^Qfc and hence the divergence of the 
k component of the coordinate qk between two parallel faces 1 and 2 at 
distance (5qfc from each other is given by 

div (J)^ = (gk,2 - qk,i) .6c[r^k/SV (5) 

where (5q^fc are the orthogonal surface area increments perpendicular to qk- 

Lemma 1 The Liouville expression for the divergence given by Eqs. (^-|^ is 
equivalent to that derived from computing it from Eq. (0j involving the absolute 
increments Sq^. 

Proof. Relative to the mid-point (p^, q^j, the Hamiltonian at the faces would 
take values H{pk,qk) it |^ l^Qk, with the positive sign for face 2, and the 
negative for face 1 and therefore Eq.(|5|) becomes, (where qk are considered 
fixed), as 



div (J) 



fdH{p, 



'k,Qk) 



+ 



y dpk 

ldH{pk,qk) 



dH{pk,qk) 



face 2 dpk 

ldH{pk,qk) 



face 1 



2 dpkdqk 
0, leads to 



Qk.l 



2 dpkdqk qk.2 
Hence taking limits as 5qk 

div{q)k = p. r. 
opkOqk 

and an exactly similar derivation for the generalized momenta yields 



div{p)k 



dpkdqk 



(6) 



(7) 



(8) 



which accords with the standard Liouville expression for the total divergence 
div{\r) = 0, since from Eqs.(|7[]8|), div{v) = Yl {div{q)k + div{p)k) = but 

in addition, we can infer by taking limits in the regions without conserving 
energy in L-space that the individual components are given by Eqs.(|^|^ and 
here, each component may be identically zero for typical systems, as shown 
below. We note that for the Louville equation, v is a non-stochastic, continuous 
and differ entiable vector. 

Theorem 1 // the Hamiltonian H is partitioned as a linear functional of the 
coordinates q and the momenta p with no cross-terms as H = K(p) + V^(q) 
then each component of the divergence is zero where div{p)k = div{q)k = 
for all indices k. 

because V does not contain 



Proof. This follows since /^ =9l&v 

oPkOlk opk \_dqk 

any pk coordinates, and by Maxwellian reciprocity, g°' g 



same reason, ^|g^ 
following obtains: 



d 



dK 



for the K variable. From Theorem |l] the 



Corollary 1 The general solution to the problem of partitioned Hamiltonian 
coordinates according to Theorem [J is -^ = -^ = for all k. 

2.1 Discussion of the solution to Corollary 1 

The general solution to the problem of Corollary 1 is therefore q = Ci + /3 
for the position coordinates and p = Di + 7 for the momentum coordinates 
where C and D are constant vectors, and t is the time coordinate. The 
motion is linear with fixed constant forces D, where the velocity vector is 

writtenv = ( -p. I • Since -^ = — v.Vp , it follows that the phase density 

gradient is orthogonal to the velocity, and such a solution implies equiden- 
sity contours of p that are all parallel to the velocity vector, like a sheath 
enclosing a line. This result contradicts for instance the well known Gibbs 
canonical distribution assumption ^ which assumes that the canonical den- 
sity p{p,q) = C exp(—PH) conforms to the Liouville equation together with 
most other current formulations that imagine v not to be constant. If we 

consider a perfect gas with mass m, for instance, then H = Yl pl-f^in and 

fc=i 

the stationary solution to the Liouville equation yields— ^p.p = 0. For a 
canonical ensemble with fixed energy E, the velocity p must be tangential to 
p , and is therefore tracing the surface of the hypersphere, and is therefore 
never a constant. It follows that all mathematical theorems that assume the 
ergodic hypothesis is necessarily false, and in fact the ergodic hypothesis is 
false from the viewpoint of the Liouville equation with continuous and differ- 
entiable variables for the given Hamiltonians. In particular, since the gener- 
alized velocities are constant, there is no way in which the Liouville equation 
can satisfy the celebrated recurrence theorem of Poincare and Zermelo with 
its multifarious and fundamental cosmological consequences, [^ which states 
: The number of phase points, which in their movement leave a given phase 
volume g without returning into it in the course of time, will be less than any 
noticeable fraction of the total number of phase points because no trajectory 
may loop back to its initial position because of the constant velocity solution 
of systems with Hamiltonians partitioned as in Theorem 1 . 
With solutions such as q^ = Ct + (3 , one might make the identification 
Qk ~ Pk/fn ~ constant , so that only linear motion is described (constant 
momentum seems implied) which contradicts the normal statistical develop- 
ment, where the Liouville equation is taken as the basis for approximate de- 
velopments in Fokker-Planck, Boltzmann and Brownian motion equations IQ]. 
The Boltzmann equation has also been developed from the above Liouville 
equation |^] although Eu |^ fundamentally denies the connection between 
the Liouville equation and Boltzmann's, stating that he used an irreversible 
kinetic equation instead of time-reversal invariant Liouville equation. It is 
common knowledge that the internal dynamics of the Boltzmann equation 
assumes time-reversal invariance, only the macroscopic outcome seems irre- 



versible, and Bolzmann justified this by stating explicitly (under the influence 
of the recurrence theorem of Poincare ) that it would take a very long time 
for the reversal of motion to be observed. It has been shown that in fact time- 
reversal in general is a false concept from a mathematical point of view |^. In 
fact it is common knowledge that the Boltzmann equation is derivable from 
the reduced density functions of the Liouville equation Q]. Hence, it is incor- 
rect to suppose that the Liouville equation as it stands is able to provide for 
the dynamics of a system from which macroscopic properties maybe computed 
analytically without introducing mathematical discontinuity and stochastic- 
ity to the probability density functions, which would lead to the breakdown 
of the equations themselves. We also note that the quantum version ^ was 
developed on principles of analogy with the classical result. Finally, Prigogine 
has also carried out an analysis of the Liouville equation based on its analogy 
with the Schrodinger equation to derive various results; based on the elemen- 
tary observation above, the analysis is of severely restricted validity, if not 



incorrect in general [10|. 



3 Development of a stochastic equation in 
/i-space 

In view of the above difficulties we develop from first principles another type 
of equation for describing irreversible phenomena in /x space, as with the 
Boltzmann equation. We recall that T space represents the set of coordinates 
for a particular system, implying that the motion of the point represents 
the evolution of the associated system. The entire ensemble of systems each 
with the same Hamiltonian Hg which are coupled to each other to form a 
Canonical ensemble constitutes one single system called the supersystem with 
the Hamiltonian Hs which contains all the coordinates of each of the systems 
j with the same Hamiltonian Hs in T' superspace; the super-Hamiltonian is 
therefore a functional of the coordinates 

R={{Q,,},{P,,}} (9) 

where the subscript i refers to the particle index (i = 1,2, ...N) for the systems 
each with A^ particles, and j refers to the system index where (j = 1, 2, ...A'^') 
for N' systems; Q and P are the spatial and momentum coordinates respec- 
tively. In r' superspace with coordinates R there is only one point which 
evolves, with the Hamiltonian written as 

N' N' N' 

Hs = Y.Hj + Y.Y.Ht^'"^ (10) 

j=i i=i j=i 

where Hj are the individual system Hamiltonians (with form Hs )contain- 
ing the coordinates of its particles as variables whereas H^^^ '^^'^'^ contain 



the (mixed) variables belonging to different systems of the ensemble due to 
the coupling of the systems; because of the coupling, the system trajecto- 
ries in r space can cross because the velocity gradients due to the system 
Hamiltonian Hg need not equal that due to the supersystem Hamiltonian Hs, 

^•^'•(w' ^) ^ i^Q' ^"j ^°^ ^^^ ^^^^ ^^* °^ coordinates {Q,P} belong- 
ing to a particular system. The gradients give the incremental trajectory of 
a system in F space and clearly relative to an isolated system (with no in- 
teractive term) all the system Hamiltonians are the same, implying the same 
trajectory for any given point in T space which cannot obtain if interactive 
terms are present as given in Eq. ( [To| ) . Single systems of N particles where each 
particle is interacting with its iV — 1 neighbors may be derived from the above 
ensemble of N' separate systems by setting j = 1. There is then a reduction 
from r space to the so-called /i space statistics of the single system, where 
the Non-Louville density D would give the normalized number of parti- 
cles within the volume element ApAq where (p, q) are the momentum and 
position coordinates common to all the particles, but in general, when j ^ 1, 
there is a reduction from V superspace to F system space. Since trajectories 
can cross at any point in F space, the simplified Liouville equation of Eq.(y) 
or other methods which employ continuous functions with non-probabilistic 
derivatives cannot be employed directly. For an ensemble of systems inter- 
acting with one another where the coordinates have been partitioned (with 
respect to the super-Hamiltonian and the coordinates that belong to a par- 
ticular system, as in Eqn.(9) so that each system is a point in F space, then 
the non-Liouville density D'^"'^ may be expressed generally as 

dD^-L(Pl.,..PN.,Qi.,..QN-,t) ^ gpN-L gpN-L gpN-L 

dt at 5Q '^ 9P ' 

= K(P,.,..PN.,Qi.,..QN.,t) (11) 

The velocities (P, Q) cannot be derived from the Liouville equation and 
so one must be able to independently determine the meaning of K as well 



as (P, Q) : Clearly, Eq.(ll) would represent in general form an evolution 
equation. Clearly, there is no immediate relation or identity between the 
density p = D of the original Liouville Eq. (pi) and the above D^~^. The P's 
and Q's of this interacting system are the independent variables of D , and 
one cannot use Hg (pertaining to a particular particle s since each particle 
is its own system) to derive velocities directly , as is sometimes assumed in 
the Liouville description. The non-Liouville density D'^"^ is given for N* 
interacting systems, each with N particles having coordinates (P^, Q^) as 

i j 

where (f(Pij) — Pa) = 1 if f(Pii) — P^ G dP^and zero otherwise and 

(f'(Qij) - Q'J = 1 if f'(Qi.) - Q1 G dQl and zero otherwise. (12) 



For the above Eq.([l2|) the primed coordinates belong to T space and the 
unprimed to T' superspace. The mapping function f(Pjj) maps the coordinate 
Pij of r' space to the T space P'i, and dP' is the volume element of arbitrary 
size (set by a particular algorithm which depends on the physical dimensions 
of the particles under study of the momentum space, with similar definitions 
for the f mapping function which acts on the Q' coordinates). We note 
that no Dirac delta functions are immediately implied by Eqn.(|l^). The 
mapping induced by Eq.(0) does not in any obvious way allow the Q and P 
variables of F' space to be continuously transformed to the Q' and P' variables 
of F space (which are variables of D'^~^). Such precise mapping functions 
must be outlined before progress may be made with an extended form of 
the Liouville equation which may detail full irreversible behavior. We now 
present a possible approach for the simplest case of an equation describing 
irreversible behavior in /u space leading to an alternative to the Boltzmann 
equation which is founded on elementary reasoning of stochastic processes, 
which does not accord completely with : 



(a) authors who disavow the use of equations of the Liouville kind |11|, be- 
cause we are able to derive a form which resembles the Liouville equation, 
but the variables can only be defined probabilistically 

(b) the utilization of the collision integral [|l^] because of the problems al- 
ready mentioned concerning the use of time-reversibility in the collision 
integral which are mathematically incorrect 

We resort to a statistical method which can be in principle verified com- 
putationally for all quantities mentioned, but where stochastic differentials 
only are involved, and where a stochastic calculus can easily be elaborated as 
needed to any degree of accuracy. 

Theorem 2 It is possible to derive a stochastic equation having the form of 
the Liouville equation which is more general than the standard Boltzmann 
equation in fi-space without recourse to any of the reversibility assumptions 
implicit in the Boltzmann equation method. 



Proof. We use the superensemble method described in Eqs.(ll-12), with 
the special case of Eq.(^) with j=l, where the super-Hamiltonian ffp/ in F' 
space (represented by one point in this 6N dimensional space) would rep- 
resent the motion of the N particles having the generalized Hamiltonian 
Hj,'{Pi,P2, ....P3N,Qi,Q2, ■■■■Qsn) with the P momentum and Q spatial co- 
ordinates with subscripted labels referring to the particles; the potential is 
completely arbitrary in the coordinates given, except that the Hamiltonian is 
assumed conservative; an example of this system is a collection of N particles 
for F' space, whereas the // space consists of the (p, q) coordinates common 
to all the particles. The /j, space Cartesian coordinates will be px,Py,ayid pz, 
for the momentum and qx,qy and qz for the spatial components. Then for 
finite increments Aqi, Api{i = x, y, z), which are components of the respective 



vectors Aq and Ap, the unnormalized fi density-in-phase -F(p, q) represent- 
ing the number of particles within the volume element ApAq would be such 
that the following equations obtain: 

Ap^ApyApzAq^AqyAq,F{p, q) = M(Ap, Aq, (P, Q)p, , (p, q)^) 



N 



y^ S'iPx - Px,j)-5'{py-Pyj).5'{pz-Pz,j)-5'{qx-Qx,j)-S'iqy-Qy,j).S'{q2:-Qz,j) 



i=i 



Yl n ^'iPi-P^,J)^'iQi-Q^ 
j=l i=x,y,z 



where 5' {pi — Pij) = 1 if \pi — Pij\ < Apj and zero otherwise, and 

^'{li ~ Qi,j) = 1 if ki ~ Qi.jl < Agj and zero otherwise. (13) 

The set of incremental limits {Ap;, Agj} are dependent on the particle char- 
acteristics and must be chosen as fixed quantities in any numerical simulation. 
Since M in an integer function, it is very clear in general that the limit lim ^ 

does not exist (e.g. as in the Wiener process), so that only a ratio based on 
finite values of the independent variables can be defined in the stochastic cal- 
culus for finite time increments < 6t < T for some specified T which is 
system dependent. To develop the above density , we need to define various 
quantities averaged over time T. 

Determination of average quantities:Th.e subscripted variables after the verti- 
cal bar refer to the variables that are kept constant in the averaging expression. 
{-St |p,q):instead of directly differentiating (|T3|) we write instead 

§F \ {-M(p,q,P,Q,t)+M(p,q,P+P5t,Q + Q<5t,t + 5t)} 



ApAq(-,p, ^^ 



(^U; (14) 



where P = — «o^)Q = ~3p^- Eq-(14) is can be simplified using stochastic 



dQ '^ ~ SP 



averaging, which will not be pursued here. Similarly, we may define for the 
momentum coordinates 

. . /'Ji^l \ {-M(pi,q,P,Q,t)+M(pi-F(^p,,q,P,Q,t)} 

/ UPi 

I F,M \ 

t,q ) (15) 



bpi ' I 5p: 

/6M 
\ Spi 



and a similar expression for the spatial coordinates where 

8 



_ {-M(p,g^,P,Q,t)+M(p,g^ + %,P,Q,t)} 

Sqi 

However, Eqs.(15-|T^) will not be further considered here. We define the vector 
component Vi of the average velocity v of the particles contained within the // 

jv' JV' 

E Pj,i E <ii,i 

space volume element ApAq as Vj = '-^^, — for momentum and Vi = '^^, for 

space variables where the subscripts i or j take coordinate labels {x, y, z) , 

and I refers to the particle label where Pi^i = pi^i and Qi^i = qi^i whenever 

{Pi,hQi,i) G '^P'^Q) and where A^' = FApAq. Since particle numbers are 

conserved we can write 

'^^ , \ ^ ^^- 



= -V.pJ^v', (17) 

i 

where in the right hand side of Eg. ([T7|) , the bar denotes the average of all 
the particle flow vectors pv'i, where i represents a particle within the fi- 
space volume increment ApAq and the non-Liouville density D^~^is given 
by D = p/ N = F(p, q, t)/A^ where N is the total number of particles in 
the system. Eq.(E17) then gives 

^|p,q^+Vp-^V,+/5V- J]vi=0. (18) 

Each particle i within the volume element ApAq must also have some coordi- 
nate {Pj, Qj) in r' space, so where the exact i particle trajectory is concerned, 
we can write the following entities: ^ = ■^- and ^ = ■^- and hence for any 
pair of velocities , the divergence terms cancel by the Hamiltonian result in 
Eqn.(^), and so V.Vj = and thus V. ^^ Vj = and Vp. ^^ Vj = V/o.v where 
V is the average velocity of the non-liouville density function p = F(p,q, t). 
Thus, we have the remarkable stochastic equivalent of the Liouville equation: 



It 



p,q / - -Vp.v (19) 



and at equilibrium, we infer ( -j^j y = =^ — Vp.v = 0. 

Remark 1 It is expected therefore that the average velocity trajectory of the 
particles would be orthogonal to the gradient of the density at equilibrium. 
This analysis does not support the Poincare theorem because there will be an 
infinite number of trajectories with different gradients for each point in p- 
space and therefore for any arbitrary phase volume Adp, dq, no limit to a 
unique trajectory exists. 

9 



On account of the importance of the above equation of motion, we shall 
refer to it as the evolution equation for the ix-space microcanonical system 
trajectory. The validity of this equation does not depend on any regime 
of particle density, nor on any dubious assumptions regarding Newtonian 
dynamics such as demanded by the Boltzmann "Stosszahlansatz"; it there- 
fore can be viewed as an alternative to the Boltzmann equation, and may 
be directly tested by computer simulation, and it would be of urgent in- 
terest to do this for simple cases, for instance hard spheres contained in a 
fixed volume at constant energy. The obvious generalization of the above 
would be to write down a stochastic equation for the /i-space evolution of 
a system connected to others in a canonical ensemble; this evolution equa- 
tion would constitute the ^-space canonical system trajectory and lastly, and 
not so obviously, it may be possible to write down an evolution equation 
for a system exchanging both energy and matter (zero-particle conservation) 
leading to the evolution equations for the ^-space Grand- Canonical system, 
trajectory , and these evolution equations can be tested by computer simula- 
tions. Eq.(|l3|) provides a means to test the Boltzmann H- Function accord- 
ing to this formulation for the non-Liouville particle distribution function. 

Lastly, for such complex systems, it would be of interest to link the general 
equations above to variational principles that have been proposed and applied 
to irreversible thermoelectrical and other [1^, 0] non-equilibrium systems 
which pertain to both (a) the dynamical flow vectors and (b) the stationary 
field and matter distribution of the system,, by analogy with Gibbs ' stationary 
principles for thermostatic equilibrium. It will be one purpose of future re- 
search to correlate the present work to the non-linear variational principles 
for non-equilibrium processes mentioned above where (a) has the form 

f dS = (20) 

dC 
and (b) the form 

= 0. (21) 



dS 

dC 



C/,{x} 



In Eq.(20), dS represents the entropy change for a "disintegrating" or "recov- 
erable" system |]l3[ along an actual trajectory, whereas Eq.(^) refers to the 
entropy change dS of a collection of particles - the macroparticle- subjected to 
a virtual perturbation in a thermofield subjected to constraints {x} pertinent 
for the entire system (such as particle number and volume of the entire system, 
but not necesarily the macroparticle) . Another topic for further research is the 
verification of the above proposed principles by direct computation over model 
systems. 
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